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Abstract 



Event-to-event fluctuations of the spatial patterns of the final states of high-enery 
' collisions, referred to as erraticity, are studied for the data generated by a soft-interaction 

model called ECOMB. The moments C P:q do not show simple power-law dependences 
on the bin size. New measures of erraticity are proposed that generalizes the bin-size 
dependence. The method should be applied not only to the soft production data of 
NA22 and NA27 to check the dynamical content of ECOMB, but also to other collision 
. processes, such as e + e~ annihilation and heavy-ion collisions. 

it 

' : 1 Introduction 

Ph. 

CD . 

Progress in the study of multiparticle production has recently been made in two distinct 
directions among many others. One is in finding measures of event-to-event fluctuations |E| 
that can probe the production dynamics more deeply than the conventional observables, such 
as multiplicity distribution and factorial moments || . Such measures have been referred to as 
erraticity ||, which quantifies the erratic nature of the event structure. The other direction 
is in the construction of a Monte Carlo generator, called ECOMB [||, that simulates soft 
interaction in hadronic collisions, capable of reproducing the intermittency data ||. ECOMB 
stands for eikonal color mutation branching, which are the key words of a model that is based 
on the parton model rather than the string model for low px processes. In this paper we 
combine the two, using ECOMB to generate events from which we calculate the erraticity 
measures. The result should be of considerable interest, since, on the one hand, the erraticity 
analysis of the NA22 data H is currently being carried out, and, on the other, it can motivate 
the investigation and comparison of erraticities in various different collision processes, ranging 
from e + e~ annihilation to heavy-ion collisions. 

The study of erraticity originated in an attempt to understand possible chaotic behaviors 
in quark and gluon jets JTJ, since QCD is intrinsically nonlinear. In the search for a measure 
of chaos it was realized that the fluctuation of the hadronic final states of a parton jet is the 
only observable feature of the QCD process that can replace the unpredictable trajectories 
in classical nonlinear dynamics. A multiparticle final state in momentum space is a spatial 
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pattern. Once a measure is found to quantify the fluctuation of spatial patterns, the use- 
fulness of the method goes far beyond the original purpose of characterizing chaoticity in 
perturbative QCD processes. Many problems involve spatial patterns; they can range from 
phase transition in condensed matter to galactic clustering in astrophysics. Even continu- 
ous time series can be transformed by discrete mapping to spatial patterns ||. Thus the 
erraticity analysis, which is the study of the fluctuation of spatial patterns, is more general 
than the determination of chaotic behavior. Indeed, we have applied it to the study of phase 
transition in magnetic systems by use of the Ising model || , as well as to the characterization 
of heartbeat irregularities in ECG time series @ . 

Multiparticle production at low px has always eluded first-principle calculation because 
of its nonperturbative nature. Various models that simulate the process can generate the 
average quantities, but fail in getting correctly the fluctuations from the averages [0]. In 
particular, few models can fit the intermittency data ||. To our knowledge ECOMB is the 
only one that can reproduce those data [||, (apart from its predecessor ECCO ||). Since 
that model is tuned to fit the data by the adjustment of several parameter, it is necessary 
to test its predictions on some new features of the production process. Erraticity is such a 
feature. The fluctuation of final-state patterns presents a severe test of any model. 

ECOMB includes many sources of fluctuations in hadronic collisions. In the framework 
of the eikonal formalism it allows for fluctuations in impact parameter b. For any b there 
is the fluctuation of the number /i of cut Pomerons. For any /i there is the fluctuation of 
the number v of partons. For any v the color distribution along the rapidity axis can still 
fluctuate initially. During the evolution process the local subprocesses of color mutation, 
spatial contraction and expansion, branching into neutral subclusters, and hadronization into 
particles or resonances can all fluctuate. Taken together the model can generate such widely 
fluctuating events that fitting some average quantity such as (n) or dn/dy does not explore 
the full extent of its characteristics. The dependence of normalized factorial moments F q on 
the bin size 5 usually called intermittency, probes deeper, but it is nevertheless a measure 
that is averaged over all events. Erraticity is a true measure of event-to-event fluctuation. 

2 Erraticity 

There are various ways to characterize a spatial pattern. We shall use the horizontal factorial 
moments. Given the rapidity distribution of a particular event, we first convert it to a 
distribution in the cumulative variable X ||, |TJ, in terms of which the average rapidity 
distribution dn/dX is uniform in X. We then calculate from that distribution for that event 
the normalized F q 

F q =(n(n-l)---(n-q + l))/(n) q , (1) 

where (• • •) signifies (horizontal) average over all bins, and n is the multiplicity in a bin. We 
emphasize that (1) does not involve any average over events. F q does not fully describe the 
structure of an event, since at any fixed q it is insensitive to the rearrangement of the bins. 
However, it does capture some aspect of the fluctuations from bin to bin, and is adequate 
for our purpose. 
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Since F q fluctuates from event to event, one obtains a (vertical) distribution P(F q ) after 
many events. Let the vertical average of F q determined from P{F q ) be denoted by (F q ) . 
Then in terms of the normalized moments for separate events 

$ q = F q /(F q ) v , (2) 

we can define the vertical pth order moments of the normalized qth order factorial (horizontal) 
moments 



C p* = {*$) v ■ ( 3 ) 
Erraticity refers to the power law behavior of C VA [|l], [| 

C m ocM^ , (4) 

where M is the number of bins, 1/5, and the length in X space is 1. ip q (p) is referred to 
as the erraticity exponent. If the spatial pattern never changes from event to event, P{F q ) 
would be a delta function at = 1, and C VA would be 1 at all M, p, and q, resulting in 
ip q {p) = 0. The larger ip q (p) is, the more erratic is the fluctuation of the spatial patterns. 

Since ip q (p) is an increasing function of p with increasing slope, an efficient way to char- 
acterize erraticity with one number (for every q) is simply to use the slope at p = 1, i.e. 



(5) 



It is referred to as the entropy index [|l|]. Experimentally, it is easier to determine first an 
entropy-like quantity E ? directly from $ g : 

V q = (%hi%) v , (6) 

which follows from (|3|) and 

= dC v Jdp\ v=1 , (7) 
and then to determine \i q from E g using 

provided that C Pi9 has the scaling behavior (|J). In [|I|] it is found that fi q is larger for quark 
jets than for gluon jets, indicating that the branching process of the former is more chaotic, 
or, in more words, the event-to-event fluctuation is more erratic. 

If the moments C VA do not have the exact scaling behavior in M, as in (|J), but have 
similar nonlinear dependences on M, we can consider a generalized form of scaling 

C m {M) oc g{Mf^ . (9) 

If @ is approximately valid for a common g(M) for all p and q, it then follows from (|7|) that 

£ 9 (M)oc/2 ff ln(/(M) , (10) 
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where 



(11) 

p=i 



Despite the similarity between ([5]) and ((TTJ) , /2 g is distinctly different from fi q and should 
not be compared to one another unless g(M) = M. 

If ( |10D is indeed good for a range of q values, then we expect a linear dependence of E g 
on S 2 as M is varied. Let the slope of such a dependence be denoted by u q , i.e., 

., = || . (12) 

Then we have 

A<Z = A2 u q . (13) 

A variation of this scheme that makes use of an extra control parameter r in the problem 
is considered in 0. It is found there that the entropy indices determined that way are as 
effective as Lyapunov exponents in characterizing classical nonlinear dynamical systems. 



3 Scaling Behaviors 

The erraticity analysis described above involves only measurable quantities, so it can be 
directly applied to the experimental data. The NA22 data at y/s = 22 GeV are ideally 
suited for this type of analysis, since F q fluctuates widely from event to event ||. The 
nuclear collision data, such as those of NA49, can also be studied, but pr cuts should be 
made to reduce the hadron multiplicity to be analyzed, thereby enhancing the erraticity to 
be quantified. 

Here we apply the analysis to hadronic collisions generated by ECOMB. The parameters 
are tuned to fit (n), P n , dn/dy and (F q ) of the NA22 data ||. Without any further 
adjustment of the parameters in the model we calculate C Pi(? (M), which are therefore our 
predictions for hadronic collisions at 22 GeV. The results from simulating 3 x 10 4 Monte 
Carlo events are shown on the left side of Fig. 1. The lines are drawn to guide the eye. 

From the points shown, it is clear that the dependences of C PS on M in the log-log 
plots are not very linear, especially for the more reliable cases of q = 2 and 3, where the 
statistics are higher. Thus the power-law behavior in (||) is not well satisfied. Since the 
general behaviors of C VA are rather similar in shape, we can regard Ci,2 as the reference 
that carries the typical dependence on M, and examine C VA vs C2.2 when M is varied as an 
implicit variable. The results are shown on the right side of Fig. 1. We have left out the 
highest points that correspond to the smallest bin size, since they show saturation at q > 2. 
We have also left out the points corresponding to InM = 0, since the scaling behaviors do 
not extend to the biggest bin size. The straightlines are linear fits of the points shown and 
lend support to the scaling behavior 

CptCcCZjr* . (14) 
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The slopes of the fits are x(p>l)i which are shown in Fig. 2. One may regard xip-iQ) as a 
representation of the erraticity properties of the particle production data, when there is no 
strict scaling law as in (|j). 

The behavior of \{Pi <?) exhibited in Fig. 2 can be described analytically, if we fit the 
points by a quadratic formula for each q. The result is shown by the lines in Fig. 2. Evidently, 
the fits are excellent. The properties of the smooth behaviors can be further summarized by 
their derivatives at p = 1: 

(15) 

P =i 

The values of x' q are 0.834, 2.818, 5.243 and 7.847 for q = 2, ■ ■ • , 5, and are shown in Fig. 3. 
We suggest that these values of x' q be used to compare with the experimental data. 

Although C p , q (M) do not satisfy (§), we can consider the more general form (H). If the 
same function g(M) is good enough in (^[) for all p and q, then it follows from ( |T4] ) that 

x(p,q) = $(p,q)/$M . (16) 

Using ( [n] ) we then have 

~ N = ^2) X ' q ■ (17) 

It should be noted that, whereas x'q follows only from the scaling property of ([14]) , the 
determination of ^(2,2), and therefore jl q , requires the knowledge of g(M) in 
To determine g(M), we write it in the form 

In g(M) = (In M) a . (18) 

By varying a, we can find a good linear behavior of In vs lng(M), as shown by the 
dashed line in Fig. 4 for a = 1.8. The corresponding value of ip(2,2) determined by the 
slope of the straightline fit is 0.119. Using that in (|i~7|) yields a set of values of fi q , which are 
shown in Fig. 5 by the open-circle points. In particular, we have 

jj 2 = 0.099 , (19) 

a quantity that has a separate significance below. 

We remark that in checking the validity of (|9|) for values of p and q other than 2, one can 
improve the linearity of the points for each p and q by slight adjustments of the value of a. 
If there is a range of possible g(M) that depends on p and q to yield the best fits, however 
small the variations in a may be, the scheme defeats the point of defining a universal ip(p, q). 
We thus propose that the emphasis of the erraticity analysis should be placed on (|T4]), which 
is independent of g{M), and that @ is examined only for p = 2, q = 2 so that (|T7D can be 
evaluated. 

Since jx q is distinct from fi q , we cannot compare our result on jl q with the theoretical 
values of fi q found for quark and gluon jets fl|], nor with the experimental values of \i q 
determined from pp collisions at 400 GeV/c (NA27) |10|| . 



x' q = Tp x(P, q ) 
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The values of jj q can also be determined independently by use of E g (M). From the 
definition in (^) we have calculated S g as functions of InM, as shown in Fig. 6(a). Not 
surprisingly, the dependences are not linear. However, when S 9 is plotted against £2 in Fig- 
6(b), they all fall into straightlines, except for the point corresponding to the smallest bin 
for q = 5 (which we have left out for the fit). The slopes, which give u q defined in flT2"|), are 
1.0, 3.244, 6.0, and 9.101 for q — 2, • • • , 5. They are shown in Fig. 7. If we examine ( |i~0|) for 
q = 2 only, and plot X 2 vs In g(M) with a = 1.8, as in Fig. 4, we obtain a linear behavior 
with a slope 

fl 2 = 0.095 . (20) 

This value is to be compared with that in (|19|) with only 4% discrepency. Of the two 
methods of determining fi2, this latter approach is more reliable, since the derivative in p at 
p = 1 is done analytically in the definition of E g in (|^), whereas in the former approach the 
differentiation is done in ( |i~5"D using the fitted curve in Fig. 2. Substituting (|20| ) into fllBf), we 
can determine the values of jx q for q > 2 from the values of uj q in Fig. 7. The result is shown 
by the solid points in Fig. 5. Clearly, the two methods yield essentially the same result. 

Another way to check the degree of consistency of the two methods, independent of the 
details on g(M), is to examine the ratio r q = y^JoJq. The quantities in that ratio are derived 
from the straightline fits of lnC p (? vs In 6*2,2 and E 9 vs S 2 (as in Fig. 1 and Fig. 6) without 
resorting to such equation as (18). According to (13) and (17), the ratio r q should be a 
constant, independent of q. From the values of x' q an d u) q given above in connection with 
Figs. 3 and 7, we find that r q = 0.834, 0.867, 0.874, and 0.862 for q = 2, 3, 4, 5. The average 
is 0.86, so the standard deviation is at the 1-2% level. Evidently, the two methods are 
quite consistent, whatever g(M) may be. From (13) and (17), one would expect r q to be 
fa/ if) (2, 2), which according to the numbers given in Fig. 4, is 0.798. The discrepency from 
0.86 is nearly 7%. Thus the disagreement of the values of jl q in Fig. 5, though not large, 
has the same root as the disagreement between (19) and (20), namely, the necessity to use 
a specific form of g(M). Nevertheless, at the level of inaccuracy of 4%, which is comparable 
to the typical uncertainty in the experimental data, the value of fa given by either (19) or 
(20) clearly provides an effective measure of erraticity in soft production. 

4 Conclusion 

In conclusion, we recapitulate the two essential points of this paper. One is the prediction 
of ECOMB on the nature of fluctuations of the factorial moments F q from event to event. 
The other is the proposed method of summarizing the scaling behaviors of C Piq that do not 
have strict power-law dependences on the bin size. The two aspects of this paper converge 
on the new erraticity measures x(p, <?)> X q i an d /V 

It is hoped that the data from both NA22 and NA27 can be analyzed in terms of these 
measures so that the dynamics of soft interaction contained in ECOMB can be checked by 
the experiments. 

The proposed measures of erraticity are, of course, more general than the application 
made here to soft production. Event-to-event fluctuation has recently become an important 
theme in collisions of all varieties: e + e~ annihilation, leptoproduction, hadronic collisions at 
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very high energies where hard subprocesses are important, and heavy- ion collisions. What 
was lacking previously is an efficient measure of such fluctuations. The erraticity measures 
proposed in [I], [|, now generalized to x{P-> <?)? Xqi u q an d P>q are we U suited for that purpose. 
They may be redundant, if strict scaling in M is good enough to give the erraticity indices 
ip(p, q). The method of treating less-strict scaling properties proposed here may well be more 
generally applicable to the wide range of collision processes amenable to erraticity study. 
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Figure Captions 

Fig. 1 Log-log plots of C Ptq versus M on the left side and versus 6*2,2 on the right side. 
The lines on the left side are to guide the eye, while the ones on the right side are 
linear fits. 

Fig. 2 The slopes of the linear fits on the right side of Fig. 1 are plotted against p for 
various values of q. The lines are fits by quadratic formula. 

Fig. 3 The derivatives of x{p, <?) hi Fig. 2 at p — 1. 

Fig. 4 The open circles are for C 2i2 and the solid points are for S 2 . The lines are linear 
fits, whose slopes are '0(2,2) and /2 2 , respectively. 

Fig. 5 jl q determined in two different ways: Eq. (17) for the open circles and Eq. (13) for 
the solid points. 

Fig. 6 (a) H q vs In M for various q; (b) S g vs £2 with the lines being linear fits. 
Fig. 7 The slopes of the straightlines in Fig. 6(b), u q , plotted against q. 
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